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Abstract. This paper addresses the problem of decomposing a numerical semigroup into m-irreducible 
numerical semigroups. The problem originally stated in algebraic terms is translated, introducing 
the so called Kunz-coordinates, to resolve a series of several discrete optimization problems. First, 
we prove that finding a minimal m-irreducible decomposition is equivalent to solve a multiobjective 
linear integer problem. Then, we restate that problem as the problem of finding all the optimal 
solutions of a finite number of single objective integer linear problems plus a set covering problem. 
Finally, we prove that there is a suitable transformation that reduces the original problem to find an 
optimal solution of a compact integer linear problem. This result ensures a polynomial time algorithm 
for each given multiplicity m. We have implemented the different algorithms and have performed 
some computational experiments to show the efficiency of our methodology. 



1. Introduction 

The rich Hterature of discrete mathematics contains an important number of references, of theo- 
retical resuhs, that have helped in solving or advancing in the resolution of many different discrete 
optimization problems. Nowadays, it is considered a standard to cite the connections between graph 
theory, commutative algebra and optimization. For instance, it is a topic to mention the connections 
between the results by Konig and Egervary in graph theory and the algorithms by Kuhn for the assign- 
ment problem or by Edmonds for the maximum matching problem f }14t il5i i20i il9t i24j). In the same 
vein, there are well-known papers that apply algebraic tools to solve single objective integer linear 
problems (see [9], [11], [22]) or multiobjective integer linear problems (see [21 IS] IH [T2]). Moreover, 
more recently we can find a rich body of literature (see [25] and the references therein) that addresses 
general discrete optimization problems (non necessarily linear) using tools borrowed from pure and 
applied Algebra (Graver bases and the like). 

On the other hand, although less known, there have been also some applications of integer program- 
ming to solve problems of commutative algebra (see [IlEllT]). The goal of this paper is to analyze and 
solve another problem arising in commutative algebra using tools from integer programming. 

A numerical semigroup is a subset S of Z+ (here Z+ denotes the set of non-negative integers) 
closed under addition, containing zero and such that Z+yS* is finite. Numerical semigroups were 
first considered while studying the set of nonnegative solutions of Diophantine equations and their 
study is closely related to the analysis of monomial curves (see 13^). By these reasons, the theory of 
numerical semigroups has attracted a number of researchers from the algebraic community. This fact 
has motivated that some of the terminology used in Algebraic Geometry has been exported to this 
field. For instance, the multiplicity, the genus, or the embedding dimension of a numerical semigroup. 
Further details about the theory of numerical semigroups can be found in the recent monograph by 
Rosales and Garcia-Sanchez 32^. 

In recent years, the problem of decomposing numerical semigroups into irreducible ones has attracted 
the interest of the research community (see [HI [27l [29l |30l [31]). Recall that a numerical semigroup 
is irreducible if it cannot be expressed as an intersection of two numerical semigroups containing it 
properly. Furthermore, more recently a different notion of irreducibility, the m-irreducibility ([6 ) has 
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appeared and has started to be analyzed. A numerical semigroup with multiplicity m is said m- 
irreducible if it cannot be expressed as an intersection of two numerical semigroups with multiplicity 
m and containing it properly. (Recall that the multiplicity of a numerical semigroup is the smallest 
non zero element belonging to it.) The question of existence of to- irreducible decompositions has been 
proved in Nevertheless, it is still missing a methodology, different to the almost pure brute force 
enumeration, to find m-irreducible decompositions of minimal size. 

In this paper, we give a methodology to obtain such a minimal decomposition into m-irreducible 
numerical semigroups by using tools borrowed from discrete optimization. To this end, we identify 
one-to-one numerical semigroups with the integer vectors inside a rational polytope (see HB )• For 
the sake of this identification, we introduce the notion of Kunz-coordinates vector to translate the 
considered problem in the problem of finding some integer optimal solutions, with respect to appro- 
priate objective functions, in the Kunz polytope. Then, the problem of enumerating the minimal 
TO-irreducible numerical semigroups involved in the decomposition is formulated as a multiobjective 
integer program. We state that solving this problem is equivalent to enumerate the entire sets of 
optimal solutions of a finite set of single-objective integer problems. The number of integer problems 
to be solved is bounded above by m — 1, where to is the multiplicity of the semigroup to be decom- 
posed. Finally, we solve a set covering problem to ensure that the decomposition has the smallest 
number of elements. Although this approach is exact its complexity is rather high and in general one 
cannot prove that it is polynomial for any given multiplicity to. This comes from the fact that there 
are nowadays relatively few exact methods to solve general multiobjective integer and linear problems 
(see [16j ) and it is known that the complexity of solving in general this type of problems is ^^P-hard. 
To overcome this difficulty, we introduce a different machinery that identify a minimal decomposition 
by solving a compact linear integer program. This approach ensures that the problem of finding a 
minimal TO-irreducible decomposition is polynomially solvable. 

In Section [2] we recall the main definitions and results needed for this paper to be selfcontained. 
Section [3] is devoted to translate the problem of finding numerical semigroups of a given multiplicity 
into the problem of detecting integer points inside a rational polytope, introducing the notion of Kunz- 
coordinates vector. We write in Section [4] the conditions, in terms of the Kunz-coordinates vector fo 
a numerical semigroup to be a TO-irreducible oversemigroup. Section [5] is devoted to formulate the 
problem of decomposing and minimally decomposing into TO-irreducible numerical semigroups as a 
mathematical programming problem. We give an exact and a heuristic approach for computing such a 
minimal decomposition based on solving some integer programming problems. In Section [6] we present 
a compact model to compute, by solving only one integer programming problem, a minimal decomposi- 
tion of a numerical semigroup into TO-irreducible numerical semigroups. There, we also prove that this 
problem is polynomially solvable. Finally, in Section [7] we show some computational tests performed 
to check the efficiency of the presented algorithms with respect to the current implementation in GAP 

m- 

2. Preliminaries 

For the sake of readability, in this section we recall the main results about numerical semigroups 
needed so that the paper is selfcontained. 

Let S" be a numerical semigroup. We say that {ni, . . . , rip} is a system of generators of if S" = 

p 

riiXi : Xi E = 1, . . . ,p}. We denote S = (ni, . . . , Up) if {tii, . . . , rip} is a system of generators 

2=1 

of S*. 

The least positive integer belonging to S is denoted by m(S'), and is called the multiplicity of S 
{ui{S) =uiin{S\{0})). 

Two important notions of irreducibility are extensively used through this paper. They are the 
following: 

Definition 1 (Irreducibility and m-irreducibility) . 
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• A numerical semigroup is irreducible if it cannot be expressed as an intersection of two nu- 
merical semigroups containing it properly. 

• A numerical semigroup with multiplicity m is to -irreducible if it cannot be expressed as an 
intersection of two numerical semigroups with multiplicity m containing it properly. 

In [B] the authors analyze and characterize the set of TO-irreducible numerical semigroups. Note 
that, in particular, any irreducible numerical semigroup is TO-irreducible, while the converse is not 
true. One of the results in that paper is the key for the analysis done through this paper and it is 
stated as follows. 

Proposition 2 ([^). Let S be a numerical semigroup with multiplicity to. Then, there exist Si, . . . , Sk 
m-irreducible numerical semigroups such that S — Si D ■ ■ ■ Cl Sk- 

From the above result one may think of obtaining the minimal number of elements involved in the 
above intersection of TO-irreducible numerical semigroup. Formally, we describe what we understand by 
decomposing and minimally decomposing a numerical semigroup with multiplicity to into m-irreducible 
numerical semigroups. 

Definition 3 (Decomposition into TO-irreducible numerical semigroups). Let S be a numerical semi- 
group with multiplicity to. Decomposing S into m-irreducible numerical semigroups consists of finding 
a set of m-irreducible numerical semigroups Si, ... , S'r(s) such that S = Si ■ ■ ■ Ci S^(s) (This decom- 
position is always possible by Proposition^^. 

A minimal decomposition of S into m-irreducible numerical semigroups is a decomposition with 
minimum t (S) (minimal cardinality of the number of m-irreducible numerical semigroups involved in 
the decomposition). 

For a numerical semigroup S, the set of gaps of S, G{S), is the set (that is finite by definition 

of numerical semigroup). We denote by g{S) the cardinal of that set, that is usually called the genus 
of S. The Frobenius number of S, F(S'), is the largest integer not belonging to S. 

Let 5 be a numerical semigroup with multiplicity m. To decompose S into m-irreducible numerical 
semigroups, we first need to know how to identify those TO-irreducible numerical semigroups. In [6] 
it is proved that S is TO-irreducible if and only if it is maximal (w.r.t. the inclusion order) in the set 
of numerical semigroups with multiplicity to and Frobenius number F(S'). In |29j it is stated that a 

"F(S') + l" 

numerical semigroup, S , is irreducible if and only if g{S) — ^ . Next, we recall two results 

in [6] that characterize, in terms of the genus and the Frobenius number, the set of m-irreducible 
numerical semigroups. 

Proposition 4 ( 6 ). A numerical semigroup with multiplicity m, S, is m-irreducible if and only if 
one of the following conditions holds: 

(1) // F(5) = g{S) = TO - 1 then 5 = {a; e Z+ : 2: > m} U {0}. 

(2) IfY(S) e {to+ 1,...,2to- 1} andg{S) = m then S ^ {x e Z+ : x > m,x ^ F(5)} U{0}. 



(3) IfF{S) > 2m then S is an irreducible numerical semigroup. (In this case g{S) 



F{ S) + 1 
2 



Corollary 5 ([6 ). Let S be a numerical semigroup with multiplicity m. Then, S is m-irreducible if 
and only if g{S) {m — l,m, — ^ — ^— !- — }. 

For a given numerical semigroup, S, our goal is to find a set of m-irreducible numerical semigroups 
whose intersection is S. Then, we can restrict the search of these semigroups to the set of numerical 
semigroups that contain S. This set is called the set of oversemigroups of S. 

Definition 6 (Oversemigroups). Let S be a numerical semigroup with multiplicity m. The set, 0{S), 
of oversemigroups of S is 

0{S) = {S' numerical semigroup : S C S'} 
The set, Om{S), of oversemigroups of S with multiplicity m is Om{S) — {S' G 0{S) : m{S') = m}. 
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Denote by Im{S) the set of minimal m-irreducible numerical semigroups, with respect to the in- 
clusion poset, in the set Om{S). From the set X„i{S) we can obtain, a first decomposition of S into 
m-irreducible numerical semigroup, although, in general, it is not minimal (see Example 27 in [6]). 

Lemma 7. Let S he a numerical semigroup with multiplicity m and Im{S) — {Si, . . . ^ Sn}- Then 
S — Si n ■ ■ ■ n Sn is a decomposition of S into m-irreducible numerical semigroups. 

Clearly, this basic decomposition is not ensured to be minimal since it may use redundant elements. 

Remark 8. Note that if S is a numerical semigroup with multiplicity m, by Proposition^ si^) — tti—I 
if and only if S = {0, m, — >■} denotes that every integer greater than m belongs to S). Hence, this 
m-irreducible numerical semigroup only appears in its own decomposition and in no one else. 

This is due to the fact that S = {0, m, — >} is the maximal element in the set of numerical semigroups 
with multiplicity m, and then Om{S) — X„i{S) = {S*} (see [6| for further details). 

From now on, we assume that S ^ S = {0, m, — >■} since by the above remark, the decomposition of 
S is trivial. 

By Proposition m and Remark [U if 5 ^ = {0,rn, — its decomposition into m-irreducible 
numerical semigroups uses two types of numerical semigroups: those that are irreducible (g(5) = 
"F(S') -|- 1" 

^ ) and those that have genus equal the multiplicity of S. 

To refine the search of the elements in I„(S'), first, we need to introduce the notion of special gap. 

Definition 9. Let S he a numerical semigroup. The special gaps of S are the elements in the following 
set: 

SG{S) = {h E G{S) : S U {h} is a numerical semigroup} 
where G{S) is the set of gaps of S . 

We denote by SGm(5') ^ {h e SG{S) : h > m}. In [B], the authors proved that S is m-irreducible 
if and only if #SGm(5') < 1 (#A stands for the cardinality of the set A). Moreover, SGm(5') = if 
and only ifS' = {0,TO,— >■} (there are no gaps greater than m in S). 

Also, if we know the special gaps of a numerical semigroup, we can search for its decomposition by 
using the following result. 

Proposition 10 ([6]). Let S, Si, . . . , Sn be numerical semigroups with multiplicity m. S — SiD- ■ -ClSn 
if and only if SG m{S) D (G(5'i) U • • • U G(S'„)) = SG,„(S'). 

From the above proposition, even if the minimal m-irreducible numerical semigroups are known, 
Tm{S) = {5*1, . . . , S*™}, some of these elements may be discarded when looking for a minimal m- 
irreducible decomposition, by checking if there are redundant elements in the intersection SG„i(5) fl 
(G(5i)U---UG(5„)). 

Then, the key is to choose elements in X,n{S) that minimally cover the special gaps of S. To this 
for, we may solve a problem fixing each of the special gaps to be covered. Note that an upper bound 
of the number of problems to be solved is the number of special gaps of a numerical semigroup that is 
bounded above by m — 1 (see |32)V 

Lemma 11. Let S ^ {0, m, — s>} be a numerical semigroup with multiplicity m, and h G SGm{S). Then, 
there exists a minimal decomposition of S into m-irreducihle numerical semigroups, S — Si D ■ ■ ■ H Sn 
such that, either h = F{Si) for some i or h ^ Si for some i such that there exists h' £ SGm('S'i) with 
Y{Si) ^h' >h. 

Proof. By Proposition [2l there exists a minimal decomposition of S into m-irreducible numerical 
semigroup, S = SiCi- ■ -OSk- By applying Proposition[TUl this decomposition must verify that SGm(S')n 
(G(5i) U • • • U G(S'„)) SG,„(5). Each special gap h e SG„^{S) must be in G{S^) for some i = 
l,...,n. Assume that h ^ F{S^) and that for all h' e SG™(50 with h' > h, F{S^) ^ h' . Then, S'^ = S^U 
{F{Si)} is a m-irreducible numerical semigroup such that SG„i(S') fl (G(5i) U • • • G(S'-) • • • U G(S'„)) = 
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SGm{S). Then, we have obtained a different minimal decomposition. (Note that it has the same 
nmxiber of terms than the original one.) 

By repeating this procedure for each h £ SG„i(S') while it is possible, we find a minimal decompo- 
sition of S fulfilling the conditions of the lemma. □ 

3. The Kunz-coordinates vector 

The analysis done through this paper uses mathematical programming tools to solve the problem of 
decomposing a numerical semigroup into m-irreducible numerical semigroups. For the sake of trans- 
lating the problem to a discrete optimization problem, we use an alternative encoding of numerical 
semigroups different from the system of generators. We identify each numerical semigroup with multi- 
plicity m with a nonnegative integer vector with m — 1 coordinates, where m is the multiplicity of the 
semigroup. To describe this identification we first need to give the notion of Apery set of a numerical 
semigroup with respect to its multiplicity 

Definition 12. Let S be a numerical semigroup with multiplicity m. The Apery set of S with respect 
to m is the set Ap(S', to) — {s € S : s — m ^ S} . 

However we are interested in the following characterization of the Apery set (see [32]): Let 5 be a 
numerical semigroup with multiplicity m, then Ap(S', m) = {0 = wq, wi, . . . , Wm-i], where Wi is the 
least element in S congruent with i modulo to, for i — 1,...,to — 1. 

Moreover, the set Ap(S', m) completely determines 5, since S — (Ap(S', to) U {to}) (see [26]), and 
then, we can identify S with its Apery set with respect to its multiplicity. Moreover, the set Ap(5, to) 
contains, in general, more information than an arbitrary system of generators of S. For instance, 

Selmer in [33] gives the formulas, g{S) = ^ (Sil.gAp(s m) ^) ~ '^^3"'" ^"^^ F(S') = max(Ap(5, to)) — to. 
Moreover, one can test if a nonnegative integer s belongs to S by checking if Wg (mod m) ^ s. The 
notion of Apery set is also given when we consider any n € S instead of to, rewriting then the definition 
adequately (see [32] )• Moreover, the smallest Apery set is Ap(iS', to). 

We consider an slightly but useful modification of the Apery set that we call the Kunz-coordinates 
vector. 

Definition 13 (Kunz-coordinates). Let S be a numerical semigroup with multiplicity to. //Ap(S', to) = 
{wq — 0, Wi, . . . , w„i^i}, with Wi congruent with i modulo to, the Kunz-coordinates vector of S is the 
vector X e Z™^"'^ with components xt = for i = l,...,m — 1. 

We say that x G Z™^^ is a Kunz-coordinates vector (or Kunz-coordinates, for short) if there exists 
a numerical semigroup whose Kunz-coordinates vector is x. 

From the Kunz-coordinates we can reconstruct the Apery set. If x G 7/^~^ is the Kunz-coordinates 
vector of 5, Ap(S', to) = {mxi -|-i:i = l,...,TO — 1}U {0}. Consequently, S can be reconstructed from 
its Kunz-coordinates. 

The Kunz-coordinates vectors have been implicitly used in |21| and |26| to characterize numerical 
semigroups with fixed multiplicity, and used in [5] to count numerical semigroups with a given genus. 

Furthermore, if S* is a numerical semigroup with multiplicity to and x e Z™^^ are its Kunz- 
coordinates, from Selmer's formulas, it is easy to compute its genus and its Frobenius number as 
follows: 

1=1 

• F{S) = m.axi{mxi -\- i} — m. (Clearly, if the maximum is reached in the i-th component, 
fIs) = i (mod to)) 

The following result that appears in |26] allows us to manipulate numerical semigroups with multi- 
plicity TO as integer points inside a polyhedron. 

Theorem 14 (Theorem 11 in |26)). Each numerical semigroup is one-to-one identified with its Kunz- 
coordinates. 
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Furthermore, the Kunz- coordinates vectors of the set of numerical semigroups with multiplicity m 
is the set of solutions of the following system of diophantine inequalities: 

Xi ^1 for all i G {l,...,rn — 1}, 

Xi + Xj — Xi+j ^0 for all \ ^ i ^ j m — 1, i + j ^ m — 1, 

Xi + Xj — Xi+j-m ^ — ^ for all 1 ^ i ^ j ^ m — 1, i + j > 

Xi e Z_|_ for all i G {1, . . . , m — 1}. 

From Theorem [TJ] and Sclmcr formulas, we can identify all the numerical semigroups (in terms 
of their Kunz-coordinates vector) with multiplicity m, genus g and Frobenius number F with the 
solutions of this system of diophantine inequalities: 

Xi ^1 for alliG{l,...,TO — 1}, 

Xi + Xj — Xi+j 5^0 for all 1 ^ i ^ J ^ m — 1, i + J ^ m — 1, 

Xi + Xj — Xi+j-m 5^ — 1 for all 1 ^ i ^ j ^ m — 1, i + j > m, 

m — 1 

X! =g 

i=l 

F — max {mxi + i} — m, 

i 

Xi e Z+ for alliG{l,...,TO — 1}. 

From the above formulation and Corollary [Sj the set of m-irreducible numerical semigroups is 
completely determined by the solutions of the following disjunctive diophantine system of inequalities 
and equations: 

Xi'^ \ for alH G {1, . . . , m — 1}, 

Xi + Xj — Xi+j ^ for all 1 ^ i ^ ^ m — 1 , i + J ^ m — 1 , 

Xi + Xj — Xi+j-m — 1 for all 1 ^ i ^ j ^ m — 1, i + j > m, 

m—1 

y Xi {m — 1, m, max {mXi + i} — m} 

i=l 

Xi gZ_|_ for all « G {1, . . . , m — 1}. 

The following result characterizes the set of oversemigroups of a numerical semigroup, in term of 
its Kunz-coordinates vector. 

Proposition 15. Let S be a numerical semigroup with multiplicity m and x G Z"'^^ its Kunz- 
coordinates. Then, the set of Kunz-coordinates vectors of oversemigroups of S with multiplicity m 
is: 

{x' G Z™^^ : x' is a Kunz-coordinates vector and x' < x} 
where < denotes the componentwise order in Z™~^. 

Proof. Let 5' G C'm(S'), and Ap(5',m) = {0,w'i, . . . ,w'^_i}. Let Ap{S,m) — {0,wi, . . . ,Wm-i}- 
The ith element in the Apery set is characterized of being the minimum element in the semigroup 
that is congruent with i modulo m. Thus, w'i < Wi for all i = l,...,m — 1, since S C S'. Then 
/ _ w--! ^ wi-i _ for alH = 1, . . . , m — 1. Hence, x' < x. □ 

Therefore, the oversemigroups of S can be identified with the ^'undercoordinates" of its Kunz- 
coordinates. 

For the ease of presentation, we identify a numerical semigroup with multiplicity m with an integer 
vector with m — 1 coordinates, its Kunz-coordinates. All the notions previously given for numerical 
semigroups are adapted conveniently by using the following notation. If 5 is a numerical semigroup 
and X G Z™^^ is its Kunz-coordinates vector, we denote: 
• m(x) — m(5) — m (Multiplicity of x). 
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F{x) = F{S) (Frobenius number). 

G(x) ~ G(5) = {n G Z : ma;„ (^nod m) + " (mod m) > n} (Gaps of x). 
g{x) — g{S) (Genus of x). 
SGix) = SG(S') (Special Gaps of x). 
SGm{x) = SGm(S') (Special Gaps greater than m of x). 

Um{x) ~ {x' e Z™^^ : a;' is a Kunz-coordinates vector and x' < x} (Undercoordinates of x). 
Observe that if x is the Kunz-coordinates vector of S, x' G Um{x) is univocally identified with 
an element S" G 0,n{S) (Proposition [15]) . 
• Ap(x) — Ap(S', m) — {0} U {mxi + i : i = 1, . . . , m — 1} (Apery set of x). 

Note that all the above indices and sets can be computed by using only the Kunz-coordinates vector 
of the semigroup. 

As assumed above, we consider that S ^ {0, to, — >}. In terms of the Kunz-coordinates, it is equiva- 

m— 1 

lent to say that x ^ {1, . . . ,1) G Z™"^ (or ^ Xi > m). 

i=l 

By Corollary [5] we say that a Kunz-coordinates vector, x £ is TO-irreducible if g{x) € {to, to — 



1, 



F{x) + 1 
2 



}. Furthermore, we say that x is irreducible if g(a;) — ^^^^'^^ . Hence, every irreducible 

Kunz-coordinates vector in Z™~^ is TO-irreducible, but the converse is not true in general. 

We also say that a set of Kunz-coordinates vectors, D = {x^,...,x''} C Z^~^, is a decompo- 
sition of a; £ Z™^^ into m- irreducible Kunz-coordinates vectors if the semigroups associated with 
the elements in D give a decomposition into TO-irreducible numerical semigroups of the semigroup 
identified with x. Equivalently, by Proposition [TUl D is a decomposition of a; £ Z^~^ into TO- 
irreducible Kunz-coordinates vectors if x'' is an TO-irreducible Kunz-coordinates vector and SG,„(a::) = 
SG™(a;)n(G(xi)U---UG(x'=)). 

Then, a minimal decomposition x £ Z™^^ into TO-irreducible Kunz-coordinates is a decomposition 
into TO-irreducible Kunz-coordinates, D — {x^ , . . . ,a;'°} C Z™""'^, with minimum cardinality. 

We define 

1m{x) — {x' £ Um{x) : x'is TO-irreducible and ^ a TO-irreducible Kunz-coordinates 
vector x* £ Um{x) such that x* > x'} 

Xm{x) is one-to-one identified with Im{S). 



4. m-IRREDUCIBLE KUNZ-COORDINATES VECTORS 

In this section we give necessary and sufficient conditions for a undercoordinate of a Kunz-coordinates 
vector to be TO-irreducible. 

Let X £ 'Z^~^ be a Kunz-coordinates vector. By the above definition, a Kunz-coordinates vector 
x' £ Um{x) if and only if there exists y £ Z™~^ such that x' + y = x. 

By applying Theorem [T4l to x' = x — y, y G Z™^^ must verify the following inequalities: 

J/i < Xi - 1 for ah I £ {1, . . . , TO - 1}, 

Vi + Vj — Vi+j ^ Xi + Xj — Xi-i-j for all 1 ^ I ^ ^ m — 1, i + j ^ m — 1, 
Ui + Vj ~ Vi+j ^ Xi + Xj - Xi+j + 1 for aU 1 < i < j ^ m — 1, i + j > m 

Actually, if we are searching for those x' — x — y that are identified with a set of m-irreducible 
undercoordinates decomposing x, we can restrict, by Corollary [SJ to consider those with genus m. 
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Therefore, y must be a solution of the fohowing system: 



(P™(a;)) 
(1) 



yi ^ Xi — I for ah i G {1, . . . , m — 1}, 

Hi + Vj — Vi+j ^ Xi + Xj — Xi^j for aU 1 ^ i ^ ^ m — 1, i + j ^ m — 1, 
Vi + Uj — Vi+j ^ Xi + Xj — Xi+j + 1 for aU 1 ^ i ^ J ^ m — 1, j + j > m, 



m — 1 

E 

1=1 



yi e M(a;,y), 



m — 1 m — 1 

Xi — m, 
1=1 i=l 



where M(x, y) — { — m, Xi — m + 1, 



i{m(a;i - y^) + i} - ?7i + 1 



}• 



RecaU that the Kunz-coordinates vector (1, . . . , 1) e Z™ ^ is not considered because it corresponds 
to S* = {0,m, — >■} that is m-irreducible, and then, its minimal decomposition is itself (Remark [5]). 

m— 1 

Clearly, these coordinates are the unique solution of the above system when constraint ([T]) is yi = 

i=l 

m— 1 

Xi - rn. 

!=1 

In the next subsections we analyze the remaining two cases for the disjunctive constraint ([Ij. 



4.1. m-irreducible undercoordinates that are irreducible. Let x £ 



be a Kunz-coordinates 



vector. First, we want to find those m-irreducible undercoordinates of x that are also irreducible. Then, 



in system (P™(a:) I, equation ([T]) is 



(2) 



i=l 



4=1 



maxi{m(xi — yi) + i} — m + 1 



Denote now by H™(x) — {y £ 



maxi{m(xi - yi) -|- i} = m{xk - yk) + k}, and by P™(a;) = 



P(x)nH™(x) for all fc = 1, . . . , m — 1. Note that H™(x) is the hyperplane in R™^^ where the Frobenius 
number oi x — y is reached in the kth component (recall that F(x) — maxjmxi + i} — m), that is, 
F(x -y) = m{xk - yk) + k - m. 

With this assumptions, P™(a;) can be described by the following system of inequalities: 



(Pfix)) 



yi ^Xi — 1 for all i 6 {1, , 
Vi+Uj - yi+] ^Xi+ Xj - 



. , m 



1}, 

for aWl^i^j^m — 1, i+j^m — 1, 
Ui+Vj — Vi+j ^ Xi + Xj — Xi+j + 1 for all 1 ^ i ^ ^ m — 1, i + j > m. 



m— 1 m— 1 

Y.y- = E- 

--^m — 1 

V gZ+ ■ 



m{xk ~ yk) + k - m + I 



Or equivalently (using that z < \z\ < z + 1 for any z G M) by a linear system of inequalities as: 
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Vi ^ Xi ~ 1 for all i G {1, . . . , m — 1}, 

Ui + Vj - Vi+j Xj - Xi+.j for all 1 ^ i < j ^ m - 1, i + j < TO - 1, 

(P™(a:)) Ui + Uj - Ui+j ^Xi+ Xj - Xi+j + 1 for all 1 < z ^ j to - 1, i + j > m, 

ni~l m — 1 

2 Ui — myk ^ 2 — mxk ~ k + m — 2, 



i=l 1=1 



2 ^Ui- myk s^2^Xi- mxk - fc + to - 1, 



4.2. m-irreducible undercoordinates with genus m. In what follows, we describe the second 
type of TO-irreducible undercoordinates of S, those with genus to. 



Xi 



TO 



m — 1 

E 

z=l i=l 

is described by the following system of diophantine inequalities: 



Denote by HG™(a;) = {y <E 



m— 1 

r'-i : ^ 



} and P™(x) = P'"(x) n HG™(x). This set 



y« + Vj 



Ui ^ Xi ~ 1 for alHG{l,...,TO — 1}, 

t/i+j ^ Xi + Xj — Xi-^-j for all 1 ^ i ^ ^ TO — 1, i + J TO — 1, 



(3) 



2/i + Vj - J/i+j =^ 



t-i+j 



1 for all 1 ^ i ^ ^ TO — 1 , i + J > TO 



m — 1 

2^1 



m — 1 

E 

1=1 



The solutions of system (P™(a;) I are easily identified by the few possible choices for the solutions of 
equation (jSj (the integer vector x — y £ Z™~^ has positive coordinates and the sum of them must be 



to). Actually, the entire set of solutions of (P™(a;) I is 



m — 1 



{x — {1, . . . ,1) — ej : Xj ^ 2, j = 1, . . . ,TO — 1} C Z 

where e^ is the jth unit vector in Z™^^. 

Then, the set of TO-irreducible undercoordinates of a: with genus to is given by the set {(1, . . . , l)+ej 
X, ^2,j = l,...,TO-l}CZ™-i. 



5. Decomposing into m-iRREDUciBLE numerical semigroups 

In the section above we characterize the TO-irreducible undercoordinates of a Kunz-coordinates 
vector X € 7/^~^ . In what follows, we use these characterizations to find a decomposition of x into 
TO-irreducible Kunz-coordinates vectors. We first give some decompositions that are not minimal in 
general by enumerating the whole set of solutions of the systems (P™(x)) and (P™(a;)). After that 



we provide a multiobjective integer programming model to obtain the set of minimal elements in 
Imix). We prove that this model is equivalent to enumerate the entire set of optimal solutions of 
some single-objective integer programming problems. Thus, a minimal decomposition can be obtained 
from the former set of solutions by solving a set covering problem. Finally, we propose a heuristic 
methodology based on the abovementioned exact approach to obtain a (minimal) decomposition of x 
into TO-irreducible Kunz-coordinates vectors. 

As a consequence of Corollary [5] and the comments above we obtain the following result that states 
how to get a decomposition into TO-irreducible Kunz-coordinates vectors by solving several systems of 
diophantine inequalities. 
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Proposition 16. Let x e Z™~ he a Kunz- coordinates vector. Any decomposition of x into m- 
irreducible Kunz- coordinates vectors is given by some elements in the form x — y where y belongs to 
the union of the solutions of the systems P™(x), . . . , P™_]^(x) and P™(x). 

Remark 17. Note that the whole set of solutions o/P™(x), . . . , P™_]^(a;) andT'^(x) gives a decompo- 
sition into m-irreducible numerical semigroups of the semigroup S identified with x. It is the maximal 
decomposition since it has the maximum possible number of m-irreducible Kunz- coordinates involved, 
all the m-irreducible undercoordinates of x. 

In the following we give a methodology to compute minimal decompositions. The main idea is to 
choose, adequately, solutions of the systems P™(x), . . . ,P™_^(x) and P™(a;). 

The first step to select decompositions that are minimal with respect to the inclusion ordering is to 
find the minimal elements in the set of m-irreducible undercoordinates of a Kunz-coordinates vector 
X. This fact can be formulated as a multiobjective integer programming problem as stated in the 
following result. 

Theorem 18. Let x G Z™^^ be a Kunz-coordinates vector. The Kunz-coordinates vectors of the 
elements in Xm{x) are in the form x — y where y is a nondominated solution of any of the following 
multiobjective integer linear programming problems. 

(MIP^'(x)) w-min(?/i,...,y„_i) s.t. y £ P^(a:), /or fc = 1, . . . , m - 1, m. 

Proof. Let x' be an element in Xm{x). Then, x' = x — y' for some y G Z™^^. If fc = F(a;') (mod m), 
then, F{x') = mx'i^-\-k — m. Since x' is an m-irreducible undercoordinate of x with the above Frobenius 
number, either y' S P™(a;) (if F(x') > 2m) or y' e P^i^) (if ^i^') < 2m). Suppose that there is a 
nondominated solution, y, of MIP™(a;) (resp. MIP™(x)) dominating y' . Then, we can find x = x — y, 
with y nondominated solution of MIP™(x) (resp. MIP™(a;)) such that y < y' and y ^ y' . Then, x > x' 
and x' ^ X, and consequently, we have found an m-irreducible maximal Kunz-coordinates in I„i{x) 
such that X > x' and x' ^ x^ contradicting the maximality oi x' . □ 

Note that, P, the union of the nondominated solutions of MIP™(x), . . . , MIP™(x) contains Im{x), 
but it may contain nondominated solutions of MIP™ {x) that dominate some nondominated solution of 
MIP™(a;), \i k ^ j. Thus, P may contain coordinates vectors that dominate one another what would 
lead to non minimal decompositions into m-irreducible Kunz-coordinates vectors. 

The key to get minimal decompositions into m-irreducible Kunz-coordinates follows by applying 
Lemma 1111 Therefore, we need to address the question about how to compute SGm(a;). Algorithm 
[1] shows the way of computing the special gaps greater than the multiplicity of a Kunz- coordinates 
vector. This algorithm is based in the following result. There, k(n) — n mod m stands for the 
nonnegative integer remainder of dividing n by m, i.e., k(n) = n mod m. 

Theorem 19. Let x £ Z™~^ be a Kunz-coordinates vector and m < h £ N. Then, h £ SGmix) if and 
only if h ~ {x]^i^ii) — 1) -)- k(/i) and such that x^|^■f -\- xj > a;i(;(k(/i)+j) ~ lk(h).j for all j = 1, . . . ,m with 

11 zf Z 7 ^ 771 

k{h) + j ^ m and 2h > mx^2h) '^^{'^^1') ! '^''T'd where jij — < „ , . for all i,j — 1, . . . ,m^l. 

Proof. The elements in SGm(x) are those elements fulfilling the following conditions (see [6]): 

• h = Wi — m, where Wi £ Ap(a;), for some i = 1, . . . , m — 1. 

• Wi — Wj ^ Ap(x) for all Wj £ Ap(a;), Wj ^ Wi. 

• 2h> w^2h) 

By the identification of Kunz-coordinates vectors and the elements in the Apery set, the first conditions 
are translated in h = mxi -\- i — m = m{xi — 1) -\- i. The second set of conditions consist of checking 
for each j ^ 7 if Wi — Wj = mxi + i — mxj — j ^ {0} U {mxk + k : k = l,...m — 1}. Note that 
if mXi -\- i — mxj — j = mxu + k for some fc, then, k(fc) — k(7 — j), so, if Wi — Wj is an element in 
Ap(a;) the unique possible choice is 7/;k(i-j). Now, if i > j, then k(7 — j) — i — j, and the condition 
is the same as checking if mxi + i — mxj — j ^ mxi^j +7 — j, equivalently, if Xi + Xi^j ^ Xj. Since 
a; is a Kunz-coordinates vector, by Theorem 1141 Xi + Xi^j > Xj, so checking that those elements are 
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different is the same that Xi + Xi-j > Xj. Clearly, by changing indices, it is that Xi + Xj > Xi+j + jij 
(in this case i + j > m). The case when i < j is analogous but taking into account that in that case 

The third conditions is direct from the algorithm given in [6] to compute SGTO(a;). □ 

The above theorem is used to compute the set SGm(a;) for any Kunz-coordinates vector x G Z"'^^ 
as shown in Algorithm [TJ 



Algorithm 1: Computing the special gaps greater than the multiplicity of a Kunz-coordinates 

vector. 

Input : A Kunz-coordinates vector x e Z™^^. 

Compute Ml — {m{xi — 1) ~\- i : Xi + xj > Xi+j, for all j with i + j < m} and 
M2 — {m{xi — I) + i : Xi + Xj > Xi+j-m — 1, for all j with i + j > to}. 

Output: SG„j(a;) = {z e Mi n M2 : z > m and 2z > mx^2z) + k(2z)}. 



Note that the complexity of Algorithm [T] is O(m^). 

From Algorithm [1] and the Kunz-coordinates vector of a numerical semigroup, we obtain the follow- 
ing useful result. 

Proposition 20. Let x G be a Kunz-coordinates vector, y € Z™~^ and h G SGm(a;)- If x — y 

is a undercoordinate of x, then, h G G(x — y) if and only if y^h) = 0. Furthermore, F{x — y) is the 
unique element in {h G SGm(a;) : k(/i) — max{i G {1, . . . , m ^ 1} : yi — 0}}. 

Proof. Since h G SG„i(a:;), by Algorithm^ h = m(xk(/i) ^ 1) + k(/i)- 

If /i G G{x - y) then, m(xk(/i) - j/k(/i)) + k(/i) >h + l = mix^h) - 1) + k(/i) + 1, that is, yk(/i) < 
m — 1 

< 1, and then yk{/i) = because yi > for alH = 1, . . . , to — 1. 

TO 

Conversely, if yk(/i) = 0, then, m{x^fi) ~ 2/k(/i)) + k(/i) — TOXk(;i) + k(/i) > h + 1 since h is an special 
gap of x, and then, in particular, a gap of x. Thus, h G G{x — y). □ 

By Proposition [TOl for each h G SGm(a;) we are looking among our solution, y, for one that holds 
h G G(a; — y). This is equivalent, by Proposition [20l to search for those with yk(/i) = 0. Then, from all 
the minimal TO-irreducible numerical oversemigroups of S, we only need for the minimal decomposition, 
those that do not contain the special gaps of S. The following result even shrink further this search. 

Lemma 21. Let x G Z™^^ be a Kunz-coordinates vector and h G SG„i(a;). Then, every nondominated 
solution of (MIP™(a;) ), y, has y^h) — 0, and then, F{x — y) = h. Moreover, y is the solution with the 
minimum sum of its coordinates (length). 

Proof. By Algorithm [TJ h — mx^h^ + k(/i) — to. Furthermore, h G G(a;') for some x' = x — y in the 
decomposition, so m{x^h) - Vkih)) + Hh) >h + l = mx^,,) + k{h) - to + 1. Then, yk(h) < 1 - ;^ < 1, 
being then yk(h) = 0. 

Then, we have a feasible solution of MIP™(a;) with ?/k(;i) = 0. Therefore, for each nondominated 
solution, y, dominating y, i.e., y < y and y ^ y. The former implies that y^h) = 0- 
Furthermore, F(a; — y) = rii{xk — yk) + k — m = mxk + k — m = h. 

Since any feasible solution, y' , of ( |MIP^(a:)[ ) must hold J2 J2 ^ Vkl + " ™ + 

i=l i=l 

. — 1 ™ 1 — 1 

fc — TO + 



then, ^ y[ 

1=1 

this length. 



i=l 



mxk 



E 

2=1 



yi, SO y has minimum length, and no one else has 

□ 



By the above result we know that, if we fix a special gap, h, a nondominated solution of MIP™(a;) 



with minimum length can be computed by fixing the value of 2/k(/i)- Then, moving through all the 
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special gaps in SGm(a;) and fixing each one of them in MIP™(x) we can obtain at least ^SGm{x) 
nondominated solutions giving a decomposition of x into m-irreducible Kmiz-coordinates. 

Therefore, an upper bound on the number of elements in any decomposition is the number of special 
gaps greater than the multiplicity of the semigroup. Thus, for each problem P^(x) we can add the 
constraint requiring that ft, is a gap of the Kunz-coordinates vector, for each h £ SGm(a:), i.e., 2/k(/i) = 0- 
Then, for each h G SG„i(a;) and k £ {!,..., m} we need to solve the following multiobjective problem: 



(MIP"(a;,ft)) 



?;-min (yi, . . . , y^.i) 
s.t. 

y^h) = 

y e Pr(^) 



Remark 22. By Lemma \21[ it is enough to search for those m-irreducible Kunz-coordinates with 
Frobenius numbers in SGmix). If h G SGm{x), this constraint is added as ma,yii{m{xi — yi)+i} — m = h, 
or equivalently as y^h) = 0- 

Note that any solution o/ MIP'"(x, /i) is a numerical semigroup with Frobenius number congruent 
with k modulo m. Since h = k{h) (mod m), h G SGm(a;), and the solutions are minimal, if one 
solution has Frobenius smaller than h, then h is not in the set of gaps of those Kunz-coordinates. 
Then, this element irrelevant for the decomposition, since there must exist some other semigroup so 
that h belongs to it. 

Hence, we can simplify further the decomposition process considering only single-objective integer 
problems rather than multiobjective ones. The following result states this fact. 

Theorem 23. Let x be a Kunz-coordinates vector. Then, the elements in a minimal decomposition of 
X into m-irreducible Kunz-coordinates must belong to the union of the set of optimal solutions of the 
following problems: 

m — l 

min ^ yi 

i=i 



(IP™(x,/i)) 



if h > 2m or 



(IP™(x,ft)) 



s.t. 



s.t. 



Vkih) = 



m — l 



yk(h) = Xk(h) 



if h < 2m, for each h G SG„i(a;). 



Proof. Problems (MIP™(x, h) I for any k = 1, . . . , m — 1 and (MIP™(a;, h)) are multiobjective programs 
with a full dimension domination cone (see 34 ). In that case, all the solutions are supported (can be 
obtained by solving scalars problems, or equivalently, the solutions are in the facets of the convex hull 
of the integer feasible region) . In our case, when fixing the special gap h, we are only interested in 
one of those nondominated solutions since all of them has h among its gaps, so they are irrelevant for 
the minimal decomposition. 

Furthermore, since the improvement cone, apart of being complete, is the cone generated by 
ei, . . . , e^-i, and then, the solutions of (IP'"(x, h) ) and (IP™(a;, h) ) are nondominated. Actually, we 
are looking for solutions with the minimum difference of gaps with x, so minimizing yi, and then, it 



is enough for our purpose to minimize the length of y as formulated in (IP'"(a:, h) ] and ( |IP™(a;, K) ]. □ 
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- 1 — ei;(ft,) ■ Furthermore, 
k(/i) — TO, the Frobenius 



Note that if (IP™(a;, h) ) is feasible, it has a unique feasible solution y = x - 
this problem is feasible if and only if k(/i) = h — m since in that case h — 2m - 
number. 

Actually, in this case, if (IP™(a;, h) ) has also a solution, y, it must be also the solution of (IP™(a;, h) ). 
It is stated in the following theorem. 

Theorem 24. Let x £ Z™^^ he a Kunz- coordinates vector, h G SG„i(a;) and and y^ optimal 
solutions of problems (IP™(.t,/i)) and (IP™(x, respectively. Then, y^ = y^ ■ 



X — y^ and x"^ ~ x — y"^ . a;^ is an 
^ is a Kunz-coordinates vector with 



Proof. We have two m-irreducible undercoordinates of x, x 
irreducible Kunz-coordinates vector with Frobenius number h. x 
Frobenius number h and genus to. Since the irreducible Kunz-coordinates are those with maximal 
genus when fixing the Frobenius number and the maximum genus in this case is m, the genus of is 



also to, since in both problems we are minimizing the length of y. 



□ 



The following result states that when solving (IP™(a;, h) I, the optimal value is known. 
Lemma 25. Let y be an optimal solution of ( |IP'"(a:,7t) I. Then, 

r I, I 1 " 



□ 



Proof. It follows directly form the satisfaction of constraint ^ and by Lemma [H] 

Let X e Z™~^ be a Kunz-coordinates vector. Once a decomposition is chosen, to select a minimal 
decomposition, we use a set covering formulation, to choose, among the overall set of minimal to- 
irreducible undercoordinates of x, a minimal number of elements for the decomposition. 

Let SGm(a;) = {hi, . . . ,hs} and Di = {x^^ , . . . , x*"* } be the set of the maximal Kunz-coordinates 
vectors of TO-irreducible undercoordinates of x when fixing the special gap hi (optimal solutions of 
IP'"(a;, hi)), for i = 1, . . . , s. We denote by £> = Z?iU- • -UDg the set of TO-irrcduciblc Kunz-coordinates 
vectors candidates to be involved in the minimal decomposition of x. 

We consider the following set of decision variables 

1 if x^^ is selected for the minimal decomposition, 
otherwise. 

for i = \, . . . ,m — \, j = \, . . . ,ip^. 

We formulate the problem of selecting a minimal number of TO-irreducible undercoordinates vectors 
of X that decompose x into TO-irreducible Kunz-coordinates as: 



(SC™(i^)) 



s Pi 

i=i ]=i 



1] 



s.t. 



E 



>1 ,V/i e SG™(x). 



;J7mx^■'^^J+k(/l)>/^^ 



The covering constraint assures that for each special gap of x there is an element in 



,,x- 



' } such that ft, is a gap of its corresponding semigroup. Minimizing the 



overall sum we find the minimum number of Kunz-coordinates fulfilling this requirement. Note that 
when solving (SC'"(-D) I at most one element in Di is choosen for each i = 1, . . . , s. 

Finally, if a numerical semigroup S with multiplicity to is given to be decomposed into TO-irreducible 
numerical semigroups, we can, by identifying it with its Kunz-coordinates, give a procedure to compute 
such a decomposition. This process is described in Algorithm [2] In that implementation we also 
consider two trivial cases: (1) when the number of special gaps greater than the multiplicity is 1, 
being then the semigroup TO-irreducible; and (2) when the number of this special gaps is 2, where the 
decomposition is given by both solutions of the two unique integer programming problems, and no 
discarding process is needed. 
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Algorithm 2: Decomposition into m-irreducible numerical semigroups. 
Input : A numerical semigroup S with multiplicity m. 

Compute the Kunz-coordinates vector of 5: x & Z™~^. (Computing the Apery set. 
D = {}. 

Compute SGTO(a;). 
if #SG„i(a;) = 1 then 
I DmIR = {x} 
else 

for hi E SG„i(a;) do 
if hi < 2m then 
I Seti:»:=£iU{l + ek(;,)}. 
else 

for each optimal solution of {IP"^ {x , h) ) , do 
|_ Set D := D U {x - f} . 
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„sl 



'}■ 



Let D ^ {x 

Let z* be an optimal solution of (SC 
Set DmIR = {x'^ e D : z*- = 1} 
Output: DmIRNS = {({m} U {ma;- +i:i = l, 



- 1}) : x' G DmIR}. 



As a consequence of all the above comments and results we state the correctness of our approach. 

Theorem 26. Algorithm computes, exactly, a minimal decomposition into m-irreducible Kunz- 
coordinates vector of a Kunz-coordinates vector x £ Z™~^. Furthermore, the entire set of solutions 



(SC™(£')) characterizes the entire set of minimal decompositions. 



Algorithm [2] is able to compute a minimal decomposition of a Kunz-coordinates vector, x G Z™ ^ , 



by enumerating the whole set of optimal solutions of (IP™(.t, /i) ). However, this task is not easy 
since, mainly, it consists of enumerating the set of solutions of a diophantine system of inequalities, 
which is hard to compute (see [?] ) . In what follows we propose an approximate approach to obtain a 



"short" decomposition into m-irreducibles by choosing an optimal solution of (IP'"(a:, ft,) ) instead of 
enumerating all of them. One may choose any of them, but we can also slightly modify the integer 
programming model to obtain a good solution. 

We consider the following set of decision variables: 

1 if h, G G(a; - y), 
otherwise. 

for i ^ 1, . . . , n, and SGm(a;) — {hi, . . . , /i„}. 

For a fixed h G SGm{x), Wi = I represents that hi is covered by the solution x — y, and then, that 
can be discarded to obtain a minimal decomposition. 

Then, to be sure that we maximize the number of elements in the previous decomposition that can 
be discarded, we formulate the problem as: 

#SG„(an) 

max Wi 

yk{h) = 

+k(/i,) - ft, - 1 + M(l - w,) > forall ft, G SG™(x) 

where M >> 
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Observe that the big-M constraint 'rn{x^f^,^ — y^%.-^) + ^{hi) — hi~l + M{1 — Wi) > assures that 
if h ^ G{x — y) (equivalently 'rn{x^f^^ — yj^^/j.-j) +k(ft,i) < hi + 1, then, Wi = 0. Otherwise, Wi could be 
or 1, but since we are maximizing, Wi = 1. 

The optimal value of this integer problem is then the number of numerical semigroups in the 
decomposition that can be discarded with this choice. 

A pseudocode of the proposed approximated scheme for obtaining a "short" decomposition of a 
Kunz-coordinates vector x G Z™^^ into m- irreducible Kunz-coordinatcs vectors by solving (IP™(x, h) I 
is shown in Algorithm [3l 



Algorithm 3: Decomposition into m-irreducible numerical semigroups. 
Input : A numerical semigroup S with multiplicity m. 

Compute the Kunz-coordinates vector of 5: x € Z™~^. (Computing the Apery set.) 
D = {}. 

Compute SGto(x). 
if #SG„i(a;) 1 then 
I DmIR = {x} 
else 

for hi E SG„i(a;) do 
if hi < 2m then 
I SetD:=Z?U{l + Ck(,,)}. 
else 

|_ Let y be an optimal solution of (lP"^{x, h) ). Set D := D U {x — i/} . 

Let D = {x\...,x'}. 
if #SG,„(.x) = 2 then 
I DmIR = D 
else 



Select a minimal decomposition from D. Let z* be an optimal solution of (SC'"(I?)). 
Set DmIR = {x^ e D : z* = 1} 

Output: DmIRNS = {({m} U {mx[ + i : i = I, . . . ,m - I}) : x' e DmIR}. 



When running Algorithm |3] we obtain an optimal solution of the problem, and then moving through 
all the special gaps we obtain a decomposition into m-irreducible Kunz-coordinates. With the following 
example we show how algorithms [2] and [3] run for a given numerical semigroup. 

Example 27. Let S = (5,11,12,18). Its Kunz-coordinates vector is x — (2,2,3,4) and SGm(S') = 
{6,13,19}. 

First, we solve one integer problem for each special gap: 

• h — 6: Since h < 2 x 5 = 10, the integer problem to solve is PKx, 6) and then Di — {x^^ = 
(2,1,1,1)}. 

• h = 13; In this case /i > 2 x 5 = 10 and h = S (mod 5), so the integer problem in this case 
is P^x, 13). The whole set of optimal solutions is {(1, 0, 0, 3), (0, 1, 0, 3)}, so D2 = {x^^ = 
(2,l,3,l),a;22 = (1,2,3,1)}. 

• Finally, forh = 19.' clearly h> 2x5= 10 and h = 4 (mod 5), so the problem is nowT'\{x, 19). 
The set of optimal solutions is {(1, 0, 0, 0), (0, 0, 0, 1)}, and then D3 = {x^^ = (1, 2, 3, 4), = 
(2,2,2,4)}. 

The above five Kunz-coordinates vectors give a decomposition in oversemigroups of S. To obtain a 
minimal decomposition we must solve the associated set covering problem. 
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Solving SC5(I?) we obtain that zn = Z31 = 1 and all other variables are set to zero, being then the 
minimal decomposition given by x^^ and x^^ , i.e., a minimal decomposition into 5-irreducible Kunz- 
coordinates is given by {(2, 1, 1, 1), (1, 2, 3, 4)}. Translating to numerical semigroups: 

S = (5,11,7,8,9) n (5,6,12,18,24) 



When solving {lP^(x,h)), we obtain the same decomposition. 



However, the decomposition obtained with Algorithm |3] may be not minimaL The foUowing example 
illustrates this fact. 

Example 28. Let S = (12,17,18,23,26,28,33,39) be a numerical semigroup with multiplicity 12. 
Its Kunz- coordinates vector is x = (4,2,3,2,1,1,3,3,2,2,1) and SGi2(x) = {21,22,27,31,32,37}. 
Then, 6 integer problems must be solved: lV^{x,21), lY'^{x,22), IP3^(a:, 27), IPy^(a;, 31), IPg^(a;, 32) 
and IV I (x, 37). By solving these problems with Xpress-Mosel 7.0 |17) we obtain the following opti- 
mal solutions x-y € {(1,1, 1,1, 1,1, 1,1,2,1, 1) , (1,1, 1,1, 1,1, 1,1, 1,2, 1), (1,2,3, 1,1, 1,1, 1,1, 1,1), 
(2,2,1,2,1,1,3,1,1,1,1), (1,2,2,2,1,1,2,3,1,1,1), (4,2,1,2,1,1,2,2,1,2,1)} 

The translations of the above coordinates in terms of numerical semigroups are 
{(12, 13, 14, 15, 16, 17, 18, 19, 20, 22, 23, 33), (12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 34, 23), 
(12, 13, 16, 17, 18, 19, 20, 21, 22, 23), (12, 15, 17, 18, 20, 21, 22, 23, 25, 26, 28, 43), 
(12, 13, 17, 18, 21, 22, 23, 26, 27,28,31,44), (12, 15, 17, 18, 21, 23, 26, 28, 31, 32, 34, 49)} 



Now, by solving problem l \SC"''{D)\ , (12,13,14,15,16,17,18,19,20,21,34,23) is discarded. Then, 
the decomposition using our methodology is given by five 12 -irreducible numerical semigroups: 
S = (12, 13, 14, 15, 16, 17, 18, 19, 20, 22, 23, 33) n (12, 13, 16, 17, 18, 19, 20, 21, 22, 23)n 
(12, 15, 17, 18, 20, 21, 22, 23, 25, 26, 28, 43) n (12, 13, 17, 18, 21, 22, 23, 26, 27, 28, 31, 44)n 
(12, 15, 17, 18, 21, 23, 26, 28, 32, 34) 

However, this decomposition is not minimal since S — (12, 13, 16, 17, 18, 19, 20, 21, 22, 23, 26, 39) H 
(12, 15, 17, 18, 20, 21, 22, 23, 25, 26, 28, 43) n (12, 13, 17, 18, 21, 22, 23, 26, 27, 28, 31, 44)n 
(12,15,16,17,18,23,26,31,32,33,34,49) is a decomposition into m-irreducible numerical semigroups 
using a smaller number of terms. 

The situation of Example [28] is due to the fact that among the whole set of optimal solutions of 



(IP™(a;, h) ], Algorithm [3] chooses a particular one, but depending of that choice more or less elements 
can be discarded from that decomposition to obtain the minimal one. To avoid this fact, we need 
to consider a compact model that connects all the possible elements in the decomposition and that 
selects, among all of them, the smallest number of solutions to decompose a Kunz-coordinates vector. 

6. A COMPACT MODEL FOR MINIMALLY DECOMPOSING INTO m-IRREDUCIBLE KUNZ-COORDINATES 

VECTORS 

In the section above we describe an exact and a heuristic procedure to compute a minimal de- 
composition of a Kunz-coordinates vector x € Z™"^ into m-irreducible Kunz-coordinates. To obtain 
solutions by using that exact procedure we need to enumerate the solutions of a knapsack type dio- 
phantine equation included in the Kunz polytope. Once we have those solutions, a set covering problem 
must be solved to obtain a minimal decomposition. By using that model, the complete enumeration 
cannot be avoided since by choosing one solution, one may obtain non-minimal decompositions when 
solving the set covering model (see Example [28| . We present here a compact model to decompose 
any Kunz-coordinates vector, x € Z™~ , merging in a single integer programming problem all the 
subproblems considered in the previous section to ensure minimal decompositions. Moreover, this ap- 
proach will allow us to prove a polynomiality result for the problem of decomposing into m-irreducible 
numerical semigoups. 

Let SGmix) = {hi, . . . , hs}. 

We consider the following families of decision variables for the new model: 

• y\ Z+, such that x — \s a. m-irreducible undercoordinate of x with Frobenius number hi, 
for alH = 1 , . . . , s and i = l,...,m — 1. 
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• wi G {0, 1}, representing ii x — y^^ is chosen (1) or not (0) for a minimal decomposition into 
m-irrcducible coordinates of x, for alH = 1, . . . , s. 



• z 



G {0, 1}, that measures if hk is a gap oi x — y (1) or not (0), for alH, fc = 1, . . . , s. Note 



that G G(a; — ?/') if and only if = 0. 

Then, the proposed model, CIP'"(a:), is described as follows: 

(CIP™(x)) 



s.t. 

(4) 
(5) 
(6) 

(7) 

(8) 

(9) 
(10) 



niin wi 
1=1 



yi<xi-i 



yl + y, 
y\ + y 



E 

z^l 
m — 1 

E 



J/i+j ^ + Xi^j 



yi+j—m ^ ^•i + m 



HE 



)wi 



= 

(hk) 



yk(hi 

E4('->i 



forall i = 1, . . . , TO — 1 forall Z = 1, . . . , s, 
if i + J < TO, forall Z = 1, . . . , s, 
if ?' + j > TO, forall Z = 1, . . . , s, 

forall ^ = 1, . . . , s with hi > 2m, 

forall / = 1, . . . , s with hi < 2m, 

forall I = 1, . . . ,s, 
forall k = 1, . . . , s. 



(11) 


^Uh>.)>^-yKhk)-M{l-wi) 


forall l,k — 


1,.. 


. ,S 


(12) 


yi^u,)<M{l^z[^^^^)) 


forall fc = 


1,.. 


. ,s 


(13) 


< wi 


forall l,k — 


1,.. 


. ,s 


(14) 


y\ G Z+, 


forall i = l,...,TO— 1,/ = 


1,.. 


. ,s 


(15) 


G {0,1}, 


forall / ~ 


1,.. 


. ,s 


(16) 


4 e {O'l} 


forall / = 1, . . . , s, j = 1, . . 


, TO 


- 1 



for Af > max{xfc(;ij) : Z = 1, . . . , s}. 

The components of any optimal solutions, y* , of the above problem in the set {y*' : y*' ^ 0,Z 



1, . . . , s} = {y*'^, . . . , J/*'''} give a minimal decomposition of x into to- irreducible Kunz-coordinates 
vectors as {a; — y*^^ : j = 1, . . . , s}. Note also that F(a; — y*'^ ) = hi.. 

Constraints (H))-® assure that x — y' is a undercoordinate of x. ^ and © give conditions related 
to the genus and the Frobenius number of those Kunz-coordinates vectors (Corollary [SJ associated to 
the choice of y' (w' = 1). Constraint © assures that hi is a gap of x — y', and (ITU)) that there is at 
least one element in the decomposition having hi among its gaps. Constraints (ITT|) - (IT^ control that 
the variables z[ are well-defined. (fT4)) - (IT6l) are the integrality and binary constraints for the variables. 



The optimal value of ( CIP™ (x) ) gives the number of Kunz-coordinates involved in a minimal de- 



composition of X into TO-irreducible Kunz-coordinates vectors. 



The solution of (CIP"'(a;)l gives exactly a minimal decomposition of x into to- irreducible Kunz- 
coordinates (or TO-irreducible numerical semigroups). However, it is harder to solve than the problems 
in Algorithm [3] since it has much more variables (by using Algorithm |3l we need to solve at most 
TO — 1 problems with to — 1 variables and a set covering problem with at most to — 1 variables while 



(CIP™(x) ) has 2(to — 1)^ + (to — 1) integer/binary variables). In the computational experiments (see 
Section we have noticed that the solutions when running Algorithm [3] are not far from minimality 



and it is faster than running Algorithm [5] or solving (CIP™(a::) I. 
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Remark 29 (m-symmetry and m-pseudosymmetry). In [6| it is also defined the notion of m- symmetry 
and m-pseudosymmetry of a numerical semigroup with multiplicity m, extending the previous notions 
of symmetry and pseudosymmetry (see [32) j. A numerical semigroup, S, with multiplicity m is m- 
symmetric if S is m-irreducible and F{S) is odd. On the other hand, S is m-pseudosymmetric if S is 
m-irreducible and F{S) is even. 

Rosales and Branco analyzed in |27j and [28j those numerical semigroups that can be decomposed 
into symmetric numerical semigroups (in this case the semigroup is called ISY-semigroup) . Another 
interesting application of this methodology is to compute a decomposition of S into m-symmetric nu- 
merical semigroups (following the notation in [28] . S is an ISYM-semigroup) . This follows by fixing in 
(CIP™(a;) ) that the m-irreducible numerical over semigroups of S associated to even special gaps do not 
appear in the decomposition (y\ — for all i = 1, . . . ,m — 1 if I is even). Thus, the m-irreducible nu- 
merical semigroups which Frobenius numbers are each one of the odd special gaps must cover the whole 
set of gaps. If this problem is feasible, its solution gives a minimal decomposition into m-symmetric 
numerical semigroups. However, in this case we cannot ensured that it is always possible to decompose 
into m-symmetric numerical semigroups (for instance, a numerical semigroup with even Frobenius 
number is not decomposable in this way). Then, if problem (CIP™(x)) is infeasible, the semigroup 
cannot be expressed as an intersection of m-symmetric numerical semigroups. 

In addition, [2 8) analyzes the set of ISYG-semigroups (those that can be expressed as an intersection 
of symmetric semigroups with the same Frobenius number). We could introduce the notion of ISYGM- 
semigroups (those that can be expressed as an intersection of symmetric numerical semigroups with the 
same Frobenius number and multiplicity). This case can be also handled with our approach by fixing 
the Frobenius number of the semigroup in (CIP™(a;) I. 

An similar methodology can be applied to compute a decomposition into m-pseudosymmetric numer- 
ical semigroups. 



Remark 30 (Computational Complexity). Assume that m is fixed. (CIP™(a;)l has at most 2(m — 
1)^ + {m — 1) variables and then, it is solvable in polynomial time ,23; . It is also worth noting that 
the heuristic approach also has polynomial time overall complexity. Indeed, for each special gap of x, 
one integer program is solved, IP™(x, /i) if h > 2m or IP™ (cc) if h < 2m. Since the number of special 
gaps is bounded above by m — I, the complexity of this step is polynomial for fixed multiplicity, so 
polynomial. Once we have the solutions for all the special gaps, the discarding step consists of solving 
the set covering problem {SC"^{D) ) with at most m — 1 variables, so polynomial in m. 

On the other hand, , the algorithm proposed in [5] to decompose a numerical semigroup S with 
multiplicity m into m-irreducible numerical semigroups can be rewritten as follows. 

Let Qx = {V,E) be a directed graph whose set of vertices is the set of undercoordinates of x, Um{x), 
and {x^,x'^) £ E if x^ — x^ — (mod m) for some h 6 SG„i(a;). Figure{li illustrates how this graph 
is built. In that figure we denote SG„i(x) — {hi, . . . ,hk} and SGrn{x + e^f^-^ — {h'l, . . . , h'f,} . The 
algorithm looks for a set o vertices {x^, . . . with the properties that #SGm(a;*) = 1 for all i — 

1, . . . , n and that any other vertex is dominated by any of the elements in the set. Furthermore, Gx is 
a tree since it does not have circuits, /n [6] , a breath first search over this tree is proposed to find the 
desired set. Clearly, the worst case complexity of this method is exponential even for fixed multiplicity. 




Figure 1. Sketch of 
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7. Computational Experiments 



In this section we present the resuhs of some computational experiments performed to analyze the 
applicability of the proposed algorithm. Our algorithm has been implemented in XPRESS-Mosel 7.0 
|17j that allows to solve the single-objective integer problems involved in the decomposition into m- 
irreducible numerical semigroups, by using a branch-and-bound method and nesting models by calling 
the library mmjobs. The algorithms have been executed on a PC with an Intel Core 2 Quad processor 
at 2x 2.50 Ghz and 4 GB of RAM. 

The complexity of the algorithm depends of the dimension of the space (multiplicity), the size of 
the coefficients of the constraints and the number of special gaps. Then, we have randomly generated 
three different batteries of numerical semigroups with the following requirements: 

Battery I: Numerical semigroups with multiplicities ranging in [0, 25] (divided in the five subin- 
tervals (0,5], (5,10], (10,15], (15,20] and (20,25]) with generators ranging in [2,5000]. Ten 
instances for each subinterval. 
Battery II: Numerical semigroups with multiplicities ranging in [10, 2000] (divided in the seven 
subintervals (10,25], (25,50], (50,100], (100,250], (250,500], (500,1000], (1000,2000]) with 
generators ranging in [2,5000]. Five instances for each subinterval. 
Battery III: Numerical semigroups with multiplicities ranging in [25, 150] (divided in the five 
subintervals (25, 50], (50, 75], (75, 100], (100, 125] and (125,150]) with generators ranging in 
[2, 5000] and with number of special gaps greater than the multiplicity less than or equal 
to 30. Ten instances for each subinterval. 
The first battery of problem is designed to compare the three algorithms: the one implemented in 



GAP, the heuristic approach (Algorithm [S]) and the compact model (CIP™(a;) ). With the second set of 



problems, we check the efficiency of Algorithm [3] for solving large instances. Finally, with the third test 



set, we compare the difficulty of solving (CIP'"(x) ) comparing to the heuristic algorithm. (Note that 



this difficulty is mainly due to the number of special gaps since it increases the number of variables.) 
Therefore, we generate numerical semigroups with very large multiplicities but where the number of 
special gaps is bounded above by 30. 

We have used recursively the function RandomListForNS of GAP 10 until we found the list of integers 
defining the semigroup with the above requirements. The implementation done for decomposing in GAP 
into m-irreducible numerical semigroups is an adaptation of the function Decomposelntolrreducibles 
for decomposing into standard irreducible numerical semigroups. 

The results of these experiments are summarized in tables [T]-[3l In these tables, m indicates the range 



of the multiplicity, CMtime and Heurtime the average times in seconds consumed by solving (CIP™(a:) I 
and Algorithm[3l respectively, in Xpress-Mosel, GAPtime informs on the average time consumed by GAP 
for the same task, #SG is the average number of special gaps of the problems and ^m-irred is the 
average number of semigroups involved in a minimal decomposition. The column avgap is the average 
difference between the number of numerical semigroups used in the heuristic decomposition and the 



number of numerical semigroups used in the minimal decomposition computed by solving (CIP'"(a;) ) 



Note that GAP was not able to solve any of the 10 instances when the multiplicity ranges in (20, 25]. 



m 


CMtime 


Heurtime 


GAPtime 


#SG 


#m-irred 


avgap 


[0,5] 


0.001 


0.020 


0.001 


1.5 


1.5 





(5,10] 


0.003 


0.054 


2.3973 


2.7 


2.3 





(10, 15] 


0.013 


0.091 


4.1645 


4.1 


3.4 


0.1 


(15,20] 


0.053 


0.081 


523.556 


5.4 


4 





(20,25] 


0.046 


0.089 


n/a 


5.7 


4.4 


0.1 



Table 1. Results of the computational experiments for Battery I. 



We have also observed that the algorithm implemented in GAP does not ensure minimal de- 
compositions into m-irreducible numerical semigroups. For instance, consider the following case: 
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m 


Heurtime 


#SG 


#m-irred 


(25,50] 


0.242 


11.8 


7.2 


(50,100] 


1.411 


19.6 


9.6 


(100,250] 


168.272 


42.4 


25.4 


(250,500] 


1318.475 


86.2 


47.8 


(500,1000] 


1056.878 


27.2 


18.8 


(1000,2000] 


1895.058 


15.2 


9.8 



Table 2. Results of the computational experiments for Battery II. 



m 


CMtime 


Heurtime 


#SG 


#m-irred 


avgap 


(25,50] 


1.064 


0.201 


9.3 


5.8 


0.7 


(50,75] 


6.981 


0.713 


13.5 


7.1 


1.1 


(75,100] 


58.580 


1.819 


16.3 


9 


1 


(100,125] 


102.999 


3.428 


15.1 


7.1 


1.6 


(125,150] 


144.531 


5.752 


15.5 


8.3 


1.3 



Table 3. Results of the computational experiments for Battery III. 



S — (15,17,19,48,52,59,73) that decomposes in GAP into six 15-irreducible numerical semigroups 
while our methodology obtains a decomposition into five 15-irreducible numerical semigroups. The 
reason why GAP fails is closely related to the same fact that prevents to ensure, in all cases. Algorithm 
[3] to get minimal solutions 

From our computational experiments we observe that except for the instances with m G [0, 5], where 
the algorithm in GAP spends almost the same time to compute the decompositions, our methodology 
solves the problems faster than GAP. Actually, in this battery solving the problem (CIP™(a::)) is the 
best way to compute such a decomposition. This is due to the minimum computational time consumed 
by Xpress-Mosel to load the problems involved in Algorithmic] 

Both, the exact algorithm based on solving (CIP™(a::)) and the heuristic approach are able to 
compute, in reasonable CPU times, minimal decompositions into m-irreducible numerical semigroups 
for multiplicities up to 150 while the procedure implemented in GAP is not able to solve problems with 
multiplicities ranging even in (20, 25]. Furthermore, although the default branch-and-bound algorithm 
is not able to solve (CIP'"(a;)l for larger multiplicities, the heuristic approach solves problems with 
multiplicities up to m = 2000. 

The heuristic approach finds, much faster than the exact approach, a short decomposition of a 
numerical semigroups into m-irreducible numerical semigroups. Furthermore, the heuristic approach 
reaches most of the times a minimal decomposition. For instance, in the first battery of problems, 
the heuristic value does not coincide with the exact optimal one in only two out the fifty instances. 
Moreover, the third battery of instances satisfies that in 30% of the cases the minimal decomposition 
coincides with the heuristic short decomposition, in 34% of the cases the difference is only one semi- 
group, in 30% of the cases is two semigroups, in 4% (two cases) is three and in only 2% (one instance) 
is four. 

Note that most of the computations done by using Algorithm [3] may be parallelized by solving in 
different cores each one of the problems (IP™(a;, h) ) since they are independent. This could improve the 
CPU times and sizes of the problems because more than 99% of the time consumed by this algorithm is 
to solving those problems, while just a little part of the time is spent solving the set covering problem. 

On the other hand, we have only implemented the proposed models in Xpress-Mosel, with the default 
branch-and-bound method. Larger instances could be solved by applying specific more sophisticated 
integer programming algorithms to solve each one of the problems. 
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